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ABSTRACT 

Aims. Stars twinkle because their light propagates through the atmosphere. The same phenomenon is expected at a longer time scale 
when the light of remote stars crosses an interstellar molecular cloud, but it has never been observed at optical wavelength. In a 
favorable case, the light of a background star can be subject to stochastic fluctuations on the order of a few percent at a characteristic 
time scale of a few minutes. Our ultimate aim is to discover or exclude these scintillation effects to estimate the contribution of 
molecular hydrogen to the Galactic baryonic hidden mass. This feasibility study is a pathfinder toward an observational strategy to 
search for scintillation, probing the sensitivity of future surveys and estimating the background level. 

Methods. We searched for scintillation induced by molecular gas in visible dark nebulae as well as by hypothetical halo clumpuscules 
of cool molecular hydrogen (H2 - He) during two nights. We took long series of 10s infrared exposures with the ESO-NTT telescope 
toward stellar populations located behind visible nebulae and toward the Small Magellanic Cloud (SMC). We therefore searched for 
stars exhibiting stochastic flux variations similar to what is expected from the scintillation effect. According to our simulations of 
the scintillation process, this search should allow one to detect (stochastic) transverse gradients of column density in cool Galactic 
molecular clouds of order of ~ 3 x 10~ 5 g/cm 2 /10000km. 

Results. We found one light-curve that is compatible with a strong scintillation effect through a turbulent structure characterized by 
a diffusion radius R di ff < \00km in the B68 nebula. Complementary observations are needed to clarify the status of this candidate, 
and no firm conclusion can be established from this single observation. We can also infer limits on the existence of turbulent dense 
cores (of number density n > 10 9 cm -3 ) within the dark nebulae. Because no candidate is found toward the SMC, we are also able to 
establish upper limits on the contribution of gas clumpuscules to the Galactic halo mass. 

Conclusions. The limits set by this test do not seriously constrain the known models, but we show that the short time- scale monitoring 
for a few 10 6 star x hour in the visible band with a > 4 meter telescope and a fast readout camera should allow one to quantify the 
contribution of turbulent molecular gas to the Galactic halo. The LSST (Large Synoptic Survey Telescope) is perfectly suited for this 
search. 

Key words. Cosmology: dark matter - Galaxy: disk - Galaxy: halo - Galaxy: structure - Galaxy: local interstellar matter - ISM: 
molecules 
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^d" \ 1 . Introduction scintillation of background stars would have a major impact on 
. the Galactic dark matter question. 
■ The present study was made to explore the feasibility of 
Q ; the detection of scintillation effects through nebulae and The SER project (Optical Scintillation by Extraterrestrial 
through hypothetical cool molecular hydrogen (H 2 - He) Refractors) is proposed to search for scintillation of extra- 
clouds. Considering the results of bary onic compact massive galactic sources through these Galactic - disk or halo - transpar- 
objects searches through microlensing O perand et al. 20071 ent h 2 clouds. This project should allow one to detect (stochas- 
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Wyrz ykowski et al. 2010[ lAlcock et al. 200(3 see also the review tic) transverse gradients of column density in cool Galactic 



of |Moniez 2010D, these clouds should now be seriously con- mo lecular clouds on the order of ~ 3 x 10" 5 g/cm 2 /10 000 km. 

sidered as a possible major component of the Galactic hid- The technique can also be used for the nebulae science. The 

den matter. It has been suggested that a hierar chical structure discovery of scintillation through visible nebulae should indeed 

of co ld H 2 could fill the Galact ic thick disk (PfennigeretaD open a new window to investigate their structure. 
1994;|Pfenniger & Combes 1994) or halo (IDe Paolis et al. 19951 

|De Paolis et al. 1998D, providmj a solution for the Galactic hid- xhe feasibility study described here concerns the search for 

den matter problem. This gas should form transparent "clumpus- scintillation through known dark nebulae such as B68 (also iden- 

cules" of ~ 30 AU size, with an average density of lO^cm" 3 , tified as LDN57 ) 9 c bl31 (also identified as B93 and L328), 

an average column density of 10 24 " 25 cm" 2 , and a surface filling through a nebula within the Circinus complex (hereafter called 

factor of ~ 1%. The detection of these structures thanks to the circinus nebula), and also includes a test for hidden matter 

search toward the SMC. As discussed below, we did not think 

Send offprint requests to: M. Moniez, e-mail: it very likely to discover a signal, and the main purpose of the 

moniez@lal . in2p3 . fr test was to predict the sensitivity of a future optical survey from 

* This work is based on observations made at the European Southern the measurement of the signal sensitivity in infrared and from 

Observatory, La Silla, Chile. the estimate of the variable star background level. 
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Fig. 1. Left: a 2D stochastic phase screen (gray scale) from a simulation of gas affected by Kolmogorov-type turbulence. 

Right: the illumination pattern from a point source (left) after crossing such a phase screen. The distorted wavefront produces 

structures at scales of ~ Rdiff (A) and R re f(A) on the observer's plane. 



2. The scintillation process 

2.1. Basics 

Refraction through an inhomogeneous transparent cloud (here- 
after called screen), which is described by a 2D phase delay 
function cp(x, y) in the plane transverse to the line of sight, dis- 
torts the wave-front of incident electromagnetic waves (Fig. 
[T])( Moniez 2003 ). The luminous amplitude in the observer's 
plane after propagation is described by the Huygens-Fresnel 
diffraction theory. For a point-like, monochromatic source, the 
intensity in the observer's plane is affected by interferences 
which take on the speckle aspect in the case of stochastic in- 
homogeneities. At least two distance scales characterize this 
speckle, which are related to the wavelength A, to the distance of 
the screen zo, and to the statistical characteristics of the stochas- 
tic phase delay function 0: 

- The diffusion radius Rdiff (A) of the screen, defined as 
the separation in the screen transverse plane for which the 
root mean square of the phase delay difference at wave- 
length A is 1 radian ( [Narayan 1992| ). Formally, Rdiff (A) is 
given in a way that < (0(x' + x,y' + y) — (p(x' ,y')) 2 >= 1 
where (x', /) spans the entire screen plane and (x, y) satis- 
fies ^x 2 + y 2 = Rdiff(A). The diffusion radius characterizes 
the structuration of the inhomogeneities of the cloud, which 
are related to the turbulence. As demonstrated in Appendix 
A, assuming that the cloud turbulence is isotropic and is de- 
scribed by the Kolmogorov theory up to the largest scale (the 
cloud's width L z ), Rdiff can be expressed as 



Rdiff(A) = 263 km x 
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cm 
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where cr 3 „ is the molecular number density dispersion within 
the cloud. In this expression, we assume that the average po- 
larizability of the molecules in the medium is a = 0.720 x 
10" 24 cm 3 , corresponding to a mixing of 76% of H2 and 24% 
of He by mass. 
The refraction radius 
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is the size in the observer's plane of the diffraction spot from 
a patch of Rdiff(A) in the screen's plane. 
- In addition, long scale structures of the screen can possi- 
bly induce local focusing/defocusing configurations that pro- 
duce long time- scale intensity variations. 



2.2. Expectations from simulation: intensity modulation, time 
scale 

After crossing an inhomogeneous cloud described by the 
Kolmogorov turbulence law (Fig[U left), the light from a 
monochromatic point source produces an illumination pattern on 
Earth made of speckles of a size Rdiff (A) within larger structures 
of a size R re f(A) (see Fig|2l up-left). The illumination pattern 
from a real stellar source of radius r s is the convolution of the 
point-like intensity pattern with the projected intensity profile of 
the source (projected radius Rs = r s x zo/zi) (Figl2l up-right). 
Rs is then another characteristic spatial scale that affects the il- 
lumination pattern from a stellar (not point-like) source. 

We simulated these illumination patterns that are caused by 
the diffusion of stellar source light through various turbulent me- 
dia as follows: 

- We first simulated 2D turbulent screens as stochastic phase 
delay functions ®(x,y), according to the Kolmogorov law 
Q; series of screens were generated at zo = 125 pc with 
100 km < R d iff < 350 km. 

- We then computed the expected illumination patterns at 
A = 2.162yum (central wavelength for K s filter) from dif- 
fused background point-like monochromatic sources located 
at z\ - I kpc using the Fast Fourier Transform (FFT) tech- 
nique (see Fig|2l up-left). Illumination patterns for other 
wavelengths and geometrical configurations were deduced 
by simple scaling. In particular, configurations compatible 
with the structure of the H2 clumpuscules of ( [Pfenniger et al.| 
1994) were produced for the / passband (Rdi ff (1. 25 fim) > 
17 km, corresponding to clumpuscules with L z = 30 AU and 
o~3n < n max = 10 10 cm" 3 ). 

- We derived series of patterns from diffused extended stellar 
sources (radii 0.257?© < r s < \.5R Q ) by convolution (Figj2] 
top-right). 

- Finally, we co-added the patterns obtained for the central and 
the two extreme wavelengths of the K s (AA = 0.275//m) and 
/ (AA = 0.290yum) passbands to simulate the diffusion of a 
wide-band source (Fig |2l down-right). 

The contrast of the patterns are clearly considerably affected by 
the size of the source (spatial coherence limitation), but only 
marginally by the bandwidth (temporal coherence). 



1 For the Kolmogorov turbulence, the 3D spectral density is a power 
law relation with exponent p = 11/3 (see appendix A). We also ex- 
plored power laws with different ft values and found the same general 
features. 
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Fig. 2. Simulated illumination map at A - 2.\6jim on Earth from 
a point source (up-left)- and from a KOV star (r s - 0.857?©, 
My = 5.9) at z\ = %kpc (right). The refracting cloud is 
assumed to be at zo = 160 pc with a turbulence parameter 
Rdi ff (2. 16 jim) - 150 km. The circle shows the projection of the 
stellar disk (with radius Rs = r s x zo/zi). The bottom maps 
are illuminations in the K s wide band (A centra i - 2A62fim, 
AA = 0.275/imJ. 



- i i i i i i i i i t < i i i i i t i i i i |. 




0.5 1 1.5 2 2.5 3 3.5 

x = R s /R ref 

Fig. 3. Expected intensity modulation index m sc i nt , = cry/7 of illu- 
mination patterns from simulated diffused stellar light as a func- 
tion of x-Rs IRref- The modulation indices are essentially con- 
tained between the curves represented by functions E m i n (x) and 
Fmax(x)- These functions allow one to constrain x when the con- 
straints on m SC i nt . are known. 



More details on this simulation will be published in a forth- 
coming paper (IHabibi et al. 201 II) . 

2.2.1. Modulation 

In general, the small speckle from a point-source is almost com- 
pletely smoothed after the convolution by the projected stellar 
profile, and only the large structures of the size R re fW - or larger 
- produce a significant modulation. Therefore, as a result of the 
spatial coherence limitations, the modulation of the illumination 
pattern critically depends on the angular size of the stellar source 
Qs - r s/(zo + Z\) ~ r s /zi. Our Monte-Carlo studies show that the 
intensity modulation index m SC i nt = o~i /I decreases when the ra- 
tio of the projected stellar disk Rs to the refraction scale R re /W 
increases, as shown in Fig. [3] This ratio can be expressed as 



Rs 



r s R diff (A) 



RrefW AZl 



-2.25 



1 jjm 



r s /zi 



R Q / 10 kpc 



RdiffW 



1000 km 



.(3) 



At the first order, as intuitively expected, the modulation 
index only depends on this ratio Rs /R re f(A) and not on other 
parameters of the phase screen or explicitly on A. Indeed, the 
dispersion of m SC i nt for series of different configurations gener- 
ated with the same Rs IRref ratio is compatible with the statis- 
tical dispersion of m SC i nt in series of 10 patterns generated with 
identical configurations. We empirically found that the modula- 
tion indices plotted as a function of Rs IRref are essentially con- 
tained between the curves represented by functions F min (x) = 

lle -l2R S IR ref Fmax ( X ) = 2e- 095R s/ R ref m Fig.^ 

2.2.2. Time scale 

Because the 2D illumination pattern sweeps the Earth with a 
constant speed, simulated light-curves of scintillating stars have 
been obtained by regularly sampling the 2D illumination pat- 
terns along straight lines. For a cloud moving with a transverse 
velocity Vt with respect to the line of sight, the velocity of the il- 
lumination pattern on the Earth is Vt(zo+zi)/zi ~ Vt aszo ^ z\\ 



neglecting the inner cloud evolution - as in radio astronomy 
(frozen screen approximation (Lyne & Graham-Smith 1998 )) - 
the flux variation at a given position is caused by this transla- 
tion, which induces intensity fluctuations with a characteristic 
time scale of 



trefW 



R re f(A) 



Vt 



' 5 .2 minutes 



(4) 



1 /j in 



1 kpc 



RdifM) 



1000 km 



V T 



100 km I s 



Therefore we expect the scintillation signal to be a stochastic 
fluctuation of the light-curve with a frequency spectrum peaked 
around l/t re f(A) on the order of (minutes)' 1 . 

As an example, the scintillation at A = 0.5 fim of a 
LMC-star through a Galactic H2 - He cloud located at 10 kpc 
is characterized by parameter Rs /R re f(0.5/im) ~ (r s /R Q ) x 
(Rdi ff (0.5 jim) /222km). According to Fig. [3] one expects 

Mscint. = (Till > 1% if X = Rs/Rref < 2.4 (= E^m)). This 

will be the case for LMC (or SMC) stars smaller than the Sun 
as soon as R^ ff (0.5 fim) < 530 km. If the transverse speed of the 
Galactic cloud is Vt = 200 km/ s, the characteristic scintillation 
time scale will be t re f ^ 24 minutes. 

2.3. Some specificities of the scintillation process 

In the subsections below we briefly describe the properties and 
specificities of the scintillation signal that should be used to dis- 
tinguish a population of scintillating stars from the population of 
ordinary variable objects. 



2.3.1. Chromaticity effect 

Because R re f depends on A, one expects a variation of the char- 
acteristic time scale t re f(A) between the red side of the optical 
spectrum and the blue side. This property is probably one of the 
best signatures of the scintillation, because it points to a prop- 
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agation effect, which is incompatible with any type of intrinsic 
source variability. 

2.3.2. Relation between the stellar radius and the modulation 
index 

As shown in Fig. [3) big stars scintillate less than small stars 
through the same gaseous structure. This characteristic signs the 
limitations from the spatial coherence of the source and can also 
be used to statistically distinguish the scintillating population 
from other variable stars. 

2.3.3. Location 

Because the line of sight of a scintillating star has to pass through 
a gaseous structure, we expect the probability for scintillation to 
be correlated with the foreground -visible gas column-density. 
Regarding the invisible gas, it may induce clusters of neighbor- 
ing scintillating stars among a spatially uniform stellar distri- 
bution because of foreground -undetected gas structures. This 
clustering without apparent cause is not expected from other cat- 
egories of variable stars. 

2.4. Foreground effects, background to the signal 

Conveniently, atmospheric intensity scintillation is negligible 
through a large telescope (m atm . 1% for a > 1 m diameter 
telescope (IDravins et al. 1997-981 )). Any other atmospheric ef- 
fect such as absorption variations at the minute scale (because of 
fast moving cirruses for example) should be easy to recognize as 
long as nearby stars are monitored together. Asteroseismology, 
granularity of the stellar surface, spots or eruptions produce vari- 
ations of very different amplitudes and time scales. A few types 
of rare recurrent var iable stars exh ibit emission variations at the 
minute scale ( Sterken & Jaschek 1996 ), but they could be identi- 
fied from their spectrum or type. Scintillation should also not be 
confused with absorption variations caused by the dust distribu- 
tion in the cloud; indeed, the relative column density fluctuations 
needed to produce measurable absorption variations (~ 1%) is 
higher by several orders of magnitudes than the fluctuations that 
are able to produce a significant scintillation (only a few 10" 7 
for the clumpuscules, and 10" 3 for the Bok globules within a 
domain of Rmff size). 

2.5. Expected optical depth 

Assuming a Galactic halo completely made of clumpuscules 
of mass M c = 10" 3 M o , their sky coverage (geometrical opti- 
cal depth) toward the LMC or the SM C should be on the or- 
der of 1% according to Pfennig er et al.| (1994). This calculation 
agrees with a simple estimat e based on the density o f the stan- 
dard halo model taken from (Caldwell & Coulson 1986). Here, 
we are considering only those structures that can be detected 
through scintillation. Therefore we quantify the sky coverage of 
the turbulent sub- structures that can produce this scintillation. 
We define the scintillation optical depth r^{Rdiff. max) as the 
probability for a line of sight to cross a gaseous cloud with a 
diffusion radius (at A) Rdiff (A) < Rdiff. max.', this optical depth is 
lower than (or equal to) the total sky coverage of the clumpus- 
cules, because it takes into account only those gaseous structures 
with a minimum turbulence strength; if Rdiff — > oo, all gaseous 
structures account for the optical depth and r^(oo) is the total sky 
coverage of the clumpuscules. 



3. Feasibility studies with the NTT 

As shown in Fig.[3j the search for scintillation induced by trans- 
parent Galactic molecular clouds makes it necessary to sample 
at the sub-minute scale the luminosity of LMC or SMC main se- 
quence stars with a photometric precision of - or better - than 
~ 1%. In principle, this can be achieved with a two meter class 
telescope with a high quantum efficiency detector and a short 
dead-time between exposures. To test the concept in a somewhat 
controlled situation, we also decided to search for scintillation 
induced by known gas through visible nebulae. 

We found that the only setup available for this short time- 
scale search was the ESO NTT-SOFI combination in the in- 
frared, with the additional benefit to enable the monitoring of op- 
tically obscured stars that are located behind dark nebulae. The 
drawback is that observations in infrared do not benefit from the 
maximum of the stellar emission and the 3.6 m diameter tele- 
scope is barely sufficient to achieve the required photometric 
precision. 

3.1. The targets 

Considering the likely low optical depth of the scintillation pro- 
cess, all our fields were selected to contain large numbers of tar- 
get stars. This criterion limited our search to the Galactic plane 
and to the LMC and SMC fields. For the search through visible 
nebulae (in the Galactic plane) we added the following require- 
ments: 

- Maximize the gas column density to benefit from a long 
phase delay. Data from 2MASS (ISkrutskie et al. 2 006) 
where used to select the clouds that induce the stronger 
reddening of background stars, pointing toward the thickest 
clouds. 

- We selected nebulae that are strongly structured (from visual 
inspection) and favored those with small spatial scale struc- 
tures. 

- We chose fields that contain a significant fraction of stars that 
are not behind the nebula to be used as a control sample. 

Our targets satisfy these requirements, except for B68, which 
does not match the second criterion, but was selected so we 
could benefit from the large number of published studies about 
this object. 

We decided to observe toward the nebulae with the K s filter 
to allow the monitoring of highly extincted stars (i.e. behind a 
high gas column density). 

To make a test search for transparent (hidden) gas, we se- 
lected a crowded field in the SMC (LMC was not observable 
at the time of observations). Then we used the / filter to col- 
lect the maximum light fluxes attainable with the SOFI detector. 
Figure |4] shows the four monitored (4.92 x 4.92) Arcmin 2 size 
fields in R and in the K s or / passbands; Table [T] gives their list 
and characteristics and also the main observational and analysis 
information. 

Through the cores of B68 and Circinus, gas column densities 
of ~ 10 22 atoms/cm 2 induce an average phase delay of 250 x 
2n at K s central wavelength (A = 2.16 fim). According to our 
studies, a few percent scintillation signal is expected from stars 
smaller than the Sun if relative column density fluctuations of 
only ~ 10" 3 occur within less than a few thousand kilometers 
(corresponding to Rdiff ^ 2000 km). These fluctuations - which 
are probably rare - could induce dust absorption variations of 
only ~ 10" 3 , which can be neglected. The expected time scale 
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Fig. 4. The four monitored fields, showing the structures of the nebulae and the background stellar densities. From left to right: B68, 
Circinus, cbl31, and the SMC. Up: images from the ESO-DSS2 in R. Down: our corresponding template images (in K s for B68, 
cbl31 and Circinus, in J for SMC). North is up and East is left. The circles on the B68 images show the position of our selected 
candidate. 





SMC 


B68 


cbl31 


Circinus 


a (J2000) 


00:52:41.3 


17:22:40.7 


18:16:59.4 


14:59:28.9 


6 (J2000) 


-72:49:14.3 


-23:49:47.2 


-18.01:53.2 


-63:06:10.1 


central gas density 




2.61 x 10 5 cm' 3 


1.8 x 10 5 cm~ 3 




center-to-edge density contrast 




16.5 


140. 




central column density Nl 




2.59 x 10 22 cm" 2 


5.8 x 10 22 cm" 2 




distance of nebula 




80 pc 


190 pc 


170 pc 


minor axis of nebula 




17000 AU 


24000 AU 


complex 


distance of sources 


62kpc 


~ 8 kpc 


~ 7 kpc 


~ 7 kpc 


light-curve duration night 1 (hours) 


2.18 


4.11 


1.25 


1.84 


light-curve duration night 2 


2.63 


5.07 


1.74 


2.16 


number of detected stars 


5042 


9599 


9084 


5249 


number of monitored stars 


691 


1114 


2779 


913 


magnitude of monitored stars 


J < 17.8 


K s < 17.1 


K s < 17.1 


K s < 17.1 


fraction of stars behind dust 


0% 


46% 


64% 




mean number of measurements/star 


980 


2013 


629 


888 



Table 1. Observations and data reduction results. The data on the nebulae are taken from Hotzel et al. (2002) for B68 and from 
\Bacma nn et al. (2000) for cbl31. The typical distance of the sources in the Galactic plane are taken from \Geo rgelin et al. (1994) 
and Russei fei al.\ ( 1 998 ). 



of the scintillation would be t re f > 5 minutes, assuming V T ~ - corresponding to the maximum stellar emission - would be 
20 km I s. better. 



3.2. The observations 

During two nights of June 2006 we took a total of 4749 consec- 
utive exposures of T exp = 10 s, with the infra-red 1024 x 1024 
pixel SOFI detector in K s (A = 2.16pm) toward B68, Circinus, 
cbl31, and in / (A = 1.25 fim) toward the SMC (see Tabled] 
for details). We recall that for a dedicated search for transpar- 
ent hidden matter, measuring visible light (B, V, R or G filters) 



4. Data reduction 

4.1. Photometric reduction 



With the EROS software JAnsari 1996b we produced the light- 
curves 0(0 of a few thousands of stars for each target. The ref- 
erence catalogs of the monitored stars where established from 
templates obtained through the standard GASGANO procedure 
(ESO 2001 ) by co-adding 10 exposures of 60s each for the B68, 
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Circinus, and cbl31 fields (13 exposures for the SMC field). 
The photometric measurements were all aligned with respect to 
these templates. We experimented with several techniques for 
the photometric reduction (aperture photometry and Gaussian 
Point Spread Function (PSF) fitting). We found that the most 
precise photometric technique, which provides the smallest av- 
erage point-to-point variation, is the Gaussian PSF fitting. We 
observed that the Gaussian fit quality of the PSF decreases when 
the seeing (i.e. the width of the PSF) is small. This effect induces 
some correlation between the seeing and the estimated flux, and 
we systematically corrected the flux for this effect, according to 
a procedure described in Tisserand (2004), Figure \5\ shows the 
dispersion of the measurements along the light-curves. This dis- 
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Fig. 5. Dispersion of the photometric measurements along the 
light-curves as a function of the mean K s magnitude (up) and 
J magnitude (down). Each dot corresponds to one light-curve. 
In the upper panel, the small black dots correspond to control 
stars that are not behind the gas; the big blue dots correspond to 
stars located behind the gas; the big star marker indicates our 
selected candidate. 



persion depends on the photometric precision and on the - pos- 
sible but rare - intrinsic stellar apparent variability. We checked 
that it is not affected by the dust by considering separately a 
control sample of stars (Fig. [7]) that are apart from the nebula 
(Fig. upper panel). We interpret the lower envelope of Fig. \5\ 
as the best photometric precision achieved on (stable) stars at a 
given magnitude. The outliers that are far above this envelope 
can be caused by a degradation of the photometric precision in a 
crowded environment, by the parasitic effect of bad pixels, or by 



a real variation of the incoming flux. The best photometric pre- 
cision does not simply result from the Poissonian fluctuations of 
the numbers of photoelectrons N ye . We discuss the main sources 
of the precision limitation in Appendix B. 



4.2. Calibration 

The photometric calibration was done using the st ars from the 
2MASS catalog in our fields (Skrutski e et al. 2006b . These stars 
were found only in the SMC and B68 fields. Because all our 
fields were observed during the same nights with stable atmo- 
spheric conditions, we extrapolated the calibration in the K s 
band from B68 to cbl3 1 and Circinus, after checking on series of 
airmass-distributed images that the airmass differences between 
the reference images (which are smaller that 0.1) could not in- 
duce flux variations higher than 0.1 magnitude. 



5. Analysis 

5.1. Filtering 

We removed the lowest quality images, stars, and mearurements 
by requiring the following criteria. 

- Fiducial cuts: We do not measure the luminosity of objects 
closer than 100 pixels from the limits of the detector to avoid 
visible defects on its borderlines. The effective size of the 
monitored fields is then 4.44 x 4.44 Arcmin 2 . 

- We remove the images with poor reconstructed seeing by 
requiring s < s + 2cr s , where s and cr s are the mean and r.m.s 
of the seeing distributions (see Fig. [6]). 

Bok Globule 68 cb131 



4500 
4000 
3500 
3000 
2500 
2000 
1500 



i 1 


m = 0.91 Arcsec. 
s = 0.13 Arcsec. . 







r 1 ^ 



n = 1.18 Arcsec. 
s =0.19 Arcsec. 



Toward the SMC 




4500 
4000 
3500 



m = 1.33 Arcsec. 
s =0.19 Arcsec. 



nJ 



X 



\ 



Seeing (Arcsec.) 



Seeing (Arcsec.) 



Fig. 6. Seeing distributions of the images towards the four tar- 
gets. The images with larger seeing than the position of the ar- 
row are discarded. 



- We remove the images with an elongated PSF that have an 
eccentricity e > 0.55. 

- We reject the measurements with a poor PSF fit quality. This 
quality varies with the seeing, the filter, the position of the 
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object on the detector, and also with the flux of the star. We 
require the^ 2 of the fit to satisfy x 1 jd.o.f < 4. 

- Minimum flux: We keep only stars with an instrumental flux 
> 1000 ADU (corresponding to J < 17.8, K s < 17.1) 
to allow a reproducibility of the photometric measurements 
better than ~ 15% in K s and ~ 8% in /. The number of these 
stars are given in Table [T] 

- We require at least 10 good quality measurements per night 
per light-curve before searching for variabilities. 

At this stage, we selected a set of sufficiently sampled and well 
measured light-curves to perform a systematic search for vari- 
abilities. The numbers of light-curves we select for each target 
are given in Table [T] As already mentioned, we divided B68 and 
cbl31 into two subfields containing control stars (out of the neb- 
ula's field) and stars behind the nebula (Fig. [7]). Table[T]gives the 
fraction of these stars located behind visible dust, which traces 
the gas toward B68 and cbl31. The case of Circinus is different 
because we cannot clearly define a control zone. 




Fig. 7. Definition of the control regions (yellow dots) and the 
search regions (blue dots) toward B68 (left) and cbl31 (right). 



global dispersion of a light-curve. R is high as soon as there is 
a correlation between consecutive fluctuations, either because of 
incoming flux variations or systematic effects. Figure [8] shows 
the distributions of R versus the apparent magnitude K s (J for 
SMC). By selecting light-curves with R > 1.6 (in K s ) or R > 1.4 
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Circinus 
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Toward the SMC 
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Fig. 8. The R = cr^/cr^ versus K s (or J) distributions of the 
light-curves of the monitored stars. The red dots correspond 
to the most variable light-curves. The green stars indicate two 
known variable objects toward SMC and our selected candidate 
toward B68. 



5.2. Selecting the most variable light-curves 

We expect the scintillation signal to produce a stochastic fluc- 
tuation of the incoming flux with a frequency spectrum peaked 
around 1 /t re f(A), on the order of (minutes)' 1 . The point-to-point 
variations of a stellar light-curve are caused by the photomet- 
ric uncertainties (statistical and systematical) and by the possi- 
ble intrinsic incoming flux fluctuations. Note that fluctuations 
with time scales shorter than the sampling (St < ti+\ -ti ~ 15 s) 
are smoothed in our data, and their high-frequency component 
cannot be detected. To select a sub-sample that includes the in- 
trinsically most variable objects, we used a simple criterion - 
not specific to a variability type -, which allowed us to both 
check our ability to detect known variable objects and explore 
new time domains of variability. Our criterion is based on the ra- 
tio R - cr^/crint of the light-curve dispersion cr^ to the "internal" 
dispersion o- int , defined as the r.m.s of the differences between 
the flux measurements and the interpolated values from the pre- 
vious and next measurements: 



(5) 



ti - tj-\ 
ti+l - ti-\ 



where N meas . is the number of flux measurements and (p(t() is the 
measured flux at time ti (the typical ti - U-\ interval is ~ 15s). 
cri nt quantifies the point-to-point fluctuations, whereas cr^ is the 



(in /) - the red dots in Fig. [8]- we retain those with a global 
variation over the two nights that is significantly larger than the 
point-to-point variations. Because we found that R almost never 
exceeds the selection threshold in large series of simulated uni- 
form light-curves affected by Gaussian errors, we systematically 
inspected each of the high R objects. We found that almost all 
of them are artifacts (all the red dots in Fig. [8]); we identified the 
following causes, which are related to the observational condi- 
tions: 

- During the meridian transit of B 6 8, the light-curves of a se- 
ries of bright stars located in two regions showed an abrupt 
flux transition, correlated with the equilibrium change in the 
mechanics of the telescope. 

- Some stars looked temporarily brighter because of contami- 
nation from the rotating egrets of bright stars. 

- We also identified a few stars transiting near hot or dead pix- 
els, which induced a distorsion of the flux determination. 

After elimination of these artifacts, only the light-curves marked 
with a star marker in Fig. [8] remain as reliable variable objects. 
For this statistically limited test, this simple way of selecting 
variable objects was sufficient, because we were able to visually 
inspect each light-curve. We also calculated the autocorrelation 
function of the light-curves with no benefit with respect to this 
basic criterion because of the artifact pollution. For future high 
statitics observations, we plan to analyze the time power spectra 
of the light-curves, which should allow an automatic segregation 
of the artifacts. 
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5.3. Sensitivity of the analysis to the known variable objects 

Our selection using variable R = 0-^/0-^ is supposed to re- 
tain any type of variability as soon as variations occur on a 
time scale longer than our sampling time. This allowed us to 
control our sensitivity to known variable stars. This control has 
been possible only toward the SMC because it is the only field 
where we found cataloged varia ble objec ts. Indeed the CDS 
and EROS cataloges (Tisserand et al. 2007 ) contain three vari- 
able objects within our SMC-field. We were able to identify 
two cepheids, #V1562 (a = 13.1550°, £ = -72.8272°, /2000, 
/ = 14.7, periodicity 4.3882 days) and the EROS cepheid 
(a = 13.2250°,^ = -72.7951°, /2000, / = 16.03, period- 
icity 2.13581 days). Figure [9] shows the corresponding folded 
light-curves from the EROS SMC database (phase diagram) 
(ITisserand et al. 20071 IHamadache et al. 20061 IRahal 20091) , on 
which we superimpose our own observations. 



Consequently, it is a star small enough to experience observ- 
able scintillation. The modulation index of the light-curve is 
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Fig. 9. EROS folded light-curves of cepheids HV1562 (upper- 
left) and of the EROS object (13.2250°, -72.7951°) (right) in 
Beros and Reros passbands, with our NTT observations in J 
(black dots). The lower-left panels show details around our ob- 
servations. 



We are content to note that the light-curves of both objects 
were successfully selected by our analysis; they correspond to 
the objects marked by a star in Fig. [H] (SMC). Our precision 
was sufficient to clearly observe the rapidly ascending phase of 
HV1562 during night 2, as can be seen in Fig. [9] 

The third variable object in the SMC field is OGLE SMC- 
SC6 148139 (a = 13.1446°, £ = -72.8333°, /2000, B = 16.5). 
It is a detached eclipsing binary with a periodicity of 1.88508 
days. Because no eclipse occured during the data taking, this 
object was - logically - not selected as a variable by our analysis. 



5.4. Signal? 

Only one variable object remains toward B68 after removing all 
artifacts (see Fig. [10| and see its location a = 260, 6762°, 6 = 
-23.8159° (/2000) in Fig.ffl-left). The star has magnitudes K s = 
16.6 and / = 20.4; its light is absorbed by dust by A K = 0.99 
magnitude; it is a main sequence star with possible type ranging 
from AO at 9.6 kpc (r s = 2.4 R Q ) to A5 at 6. 1 kpc (r s = 1 .7 R Q ) or 
from F0 at 5.0 kpc (r s = 1.6 R Q ) to F5 at 4.0 kpc (r s = IAR Q ). 



Star #4338 
t = 91.8 min. 



♦ * 



Star #4338 
t = 176.2 min 



1 



Fig. 10. Light-curves for the two nights of observation (top) and 
images of the selected candidate toward B68 during a low- 
luminosity phase (middle) and a high- luminosity phase (bot- 
tom); North is up, East is left. 



m = 0.17, which is quite high but not incompatible with a varia- 
tion owing to scintillation. If we consider that this modulation is 
caused by a scintillation effect, then it is induced by a turbulent 
structure with Rs /R re f < 0.25, according to Fig.0 as we know 
that r s /z\ > 2.5R Q / 10 kpc, we directly conclude from expression 
([3]) that Rdiff (2.1 6yum) < 96 km; using expression (Q]), this gives 
the following constraint on the size and the density fluctuations 
of the hypothetic turbulent structure: 



o~ 3n > 1.45 x 10 9 cm 3 



17000 AC/ 



(6) 



The largest possible outer scale L z corresponds here to the minor 
axis of B68 (17000 AU), but smaller turbulent structures within 
the global system may also happen. 

A definitive conclusion on this candidate would need com- 
plementary multi-epoch and multicolor observations, as ex- 
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plained in Sect. 12.3.11 but the hypothetic turbulent structure that 
is possibly responsible for scintillation has probably moved from 
the line of sight since the time of observations, considering its 
typical size. Nevertheless, reobserving this object would allow 
one to check for any other type of variability. Considering the 
short time scale and the large amplitude fluctuations, a flaring or 
eruptive star may be suspected, but probably not a spotted star 
or an effect of astroseismology. 

An important result comes from the rarity of these fluctuat- 
ing objects: There is no significant population of variable stars 
that can mimic scintillation effects, and future searches should 
not be overwhelmed by background of fakes. 

Ideally, for future observation programs aiming for an un- 
ambiguous signature, complementary multiband observations 
should be planned shortly after the detection of a scintillation 
candidate. After detection of several candidates, one should also 
further investigate the correlation between the modulation index 
and the estimated gas column density to reinforce the scintilla- 
tion case (see Sect. l2.3.2l and l2.3.3l) . 



6. Establishing limits on the diffusion radius 

In the following sections we will establish upper limits on the 
existence of turbulent gas bubbles based on the observed light- 
curve modulations. The general technique consists to find the 
minimum diffusion radius Rdiff for each monitored star that is 
compatible with the observed modulation. 

The smallest diffusion radius compatible with observed 
stellar light-curve fluctuations. 

Let us consider a star with radius r s , placed at distance z\ behind 
a screen located at zo (with a projected radius Rs = r s x zo/zi). If 
a turbulent structure characterized by Rdiff (with the correspond- 
ing R re f given by Eq. d2])) induces scintillation of the light of this 
star, the corresponding modulation m SC i nt is within the limits 



F m in(Rs I Rref) < ^scint. < F max (Rs / 'Rref), 



(7) 



as predicted from Fig. [3] 

Below we assume that our time sampling is sufficient to take 
into account any real variation within a time scale longer than 
a minute. Then, because the observed modulation m of a stellar 
light-curve results from the photometric uncertainties and from 
the hypothetic intensity modulation m SC i nt , one can infer that 
niscint. ^ m - This inequality combined with inequality O leads 
to F min (Rs I Rref) < m. Because F min is a decreasing function, 
it follows that Rs /Rref > F~^ in (m). Using Eq. © and inverting 
function F min (x) = 0A7e~ l2x , this can be expressed as a con- 
straint on the value of Rdiff for the gas crossed by the light: 



R, 



diff 



. r>min 
> K diff 



310 km 



A 



Ipm 



r s /zi 



Ro/ 10 kpc 



In 



0.17 



m 



(8) 



Information on the source size. 

To achieve a star-by- star estimate of R™^p it appears that we 
need to know zi, the distance from the screen to each source, 
and each source size r s . But because z\ » zq, it follows that 
Z\ ~ Zo + Zi, and therefore the knowledge of the angular stellar 
radius 6 S = r s /(zo + Zi) ~ r s /zi - which can be extracted from 
the apparent magnitude and the stellar type - is sufficient for this 
estimate. 

We can extract constraints on the stellar apparent radius 6 S 
from the K s (or /) apparent magnitudes by using the following 
relation derived from the standard Stefan-Boltzman blackbody 



law: 



log 



6(Rq at lOkpc) 
K s (orJ) (V-K s orJ) 



(9) 



+ 



M V o+BCo-BC 



-2 log 



T 



where R Q My©, BC Q , T Q are the solar radius, absolute V mag- 
nitude, bolometric correction and temperature; 0(R Q at lOkpc) 
is the angular solar radius at lOkpc, (V - K s or J) tak en from 
(Johnson 1966), BC and T taken from ( |Cox (Al len)) are the 
color index, bolometric correction, and temperature respectively 
that characterize the type of the star (independently of its dis- 
tance). We use this relation to establish the connection between 
K s (or /) and the angular stellar radius 6 S for a given stellar type. 
If the distance to the star is known, then its type is directly ob- 
tained from its location in the color-magnitude diagram; if the 
distance is uncertain, spanning a given domain, we estimate an 
upper value of 6 S as follows: for each stellar type, we calculate 
the distance where the apparent magnitude of the star would be 
K s . If this distance is within the allowed domain, we estimate 
6 S from Eq. (0). If the branch of the star (main sequence or red 
giant) is known, we restrain the list of types accordingly. We 
conservatively consider the maximum value &^ ax found with this 
procedure. In the next two sections, we will use Eq. ([8]) together 
with the stellar size constraints for each monitored star to estab- 
lish distributions of R™jff toward the SMC and toward the dark 
nebulae. 



7. Limits on turbulent Galactic hidden gas toward 
the SMC 

Here, the distance of t he stars are all the same zo + z\ ~ 
Z\ = 62 kpc ( Szew czyk et al. 2009] ). Therefore the angular 
stellar radii 6 S can be (roughly) estimated from the observed 
/ magnitude. We used the EROS d ata (ITisserand et al. 20071 
IHamadache et al. 20061 IRahal 20091) to distinguish between 
main sequence and red giant branch stars from the Reros ver- 
sus (B E ros - Reros ) color-magnitude diagram; then the stel- 
lar type is obtained from the absolute magnitude Mj = J - 
5. log(62000/?c/10/?c) and 6 S is derived from Eq. ([9]) (see Fig.fTTl 
left). With these ingredients, we can estimate from Eq. ® R™%f 
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Fig. 11. (left) Angular radius 6 S and type of the SMC stars as a 
function of their apparent magnitude J ( in units of the angular 
radius of the Sun at lOkpc). 

( right) Maximum angular radius 0^ ax and corresponding type of 
the B68 stars as a function of the apparent magnitude K s . The 
blue lines belong to the main sequence stars, the red lines belong 
to the red giant stars. 
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the lowest value of Rdiff compatible with the measured modu- 
lation index m for each stellar light-curve. Figure [T2l a) shows 
the cumulative distribution of this variable R™^p that is N*(Rd), 
the number of stars whose line of sight do not cross a gaseous 
structure with Rdiff < Rd- The bimodal shape of this distribu- 
tion is caused by the prominence of the red giants in the mon- 
itored population: scintillation is expected to be less contrasted 
for red giant stars because of their large radius; therefore, struc- 
tures with Rdiff k, 250 km cannot induce detectable modulation 
in the red giant light-curves, and the R^lff vames are lower than 
250 km. 

The distribution vanishes beyond Rd ~ 800 km because our 
limited resolution on the main sequence stars (a few %) pre- 
vents us from detecting any scintillation of gaseous structures 
with Rdiff > 800 km. The possible Rdiff domain for th e hidden 
gas clumpuscules expected from the model of (Pfenniger et al. 
1994; [Pfenniger & Combes 1994 ) and their maximum contribu- 
tion to the optical depth are indicated by the gray zone; the 
minimum Rdiff (~ 17 km at X - \.25pm) for these objects 
is estimated from Eq. (Q]) assuming the clumpuscule's outer 
scale is 30 AU and considering the maximum possible value 
of the density fluctuation, mathematically limited by the max- 
imum density (o~™ x < n max = 10 1Q cm~ 3 )([Pfenniger et al.| 1994; 



Pfenniger & Combes 1994) . 

From the distribution of N*(Rd), one can infer limits on the 
scintillation optical depth r^s^mifid) as it is defined in Sect. 12.51 
Indeed N*(Rd) is the number of lines of sights (l.o.s.) that do 
not cross structures with Rdiff < Rd- Defining N 'behind as the to- 
tal number of monitored l.o.s. through the nebula (all monitored 
stars in the present case because we are searching for invisible 
gas), the upper limit on the optical depth 7*1.25^ is the up- 
per limit on the ratio p = NbM if N ^ . The 95% statistical upper 

1 1 A 'behind 1-1 

limit on p is given by the classical confidence interval i.e. 0: 
- if N*(Rd) and Nbehind ~ N*(Rd) are both larger than 4, then 



T X {R d )<p + 1.643 



Pd-P) 



N> behind 

- if N*(R d ) < 4, then 

Nbehind ~ Ngs%(Rd) 



T A (Rd) < 



N, behind 



(10) 



(ii) 



where Ngs%(Rd) is the 95%C.L. Poissonian lower limit on 
- if ^behind ~ N*(R d ) < 4, then 



rx{R d ) < 



[Nbehind _ M.(/?rf)]95% 

N behind 



(12) 



where [TV 'behind ~N*(Rd)]95% is the 95%C.L. Poissonian upper 
limit on N be hind ~ N*(R d ). 

As an example, the upper limit on the optical depth for structures 
with Rdiff < 400 km is obtained from the value of Af*(400 km) ~ 
100; as Nbehind = 691, it comes p ~ 0.86. For this example, the 
95% CL upper limit on t\ 25^ (400 km) is found to be 0.88 from 
Eq. flg. 

The expected optical depth is proportional to the total mass 
of gas. In Pfenniger et al.| (1994) the clumpuscules are expected 



2 the first formula gives the upper limit of the 90% CL interval for the 
p value (see classical texbooks like Ventsel T973l ); therefore the proba- 
bility that the true p value is higher than the upper limit of this interval 
is 5%. 




500 1000 

R d (km) 

Fig. 12. (a) N*(Rd), the number of SMC stars (or lines of sight) 
with no turbulent structure of Rdiff '(!• 25pm) < Rd along the 
line of sight. The gray band shows the allowed Rdiff (1. 25pm) 
region for clumpuscules (see text), (b) The 95% CL maximum 
optical depth of structures with Rdiff < Rd toward the SMC. The 
right scale gives the maximum contribution of structures with 
Rdiff (1. 25pm) < Rd to the Galactic halo (infraction); the gray 
zone gives the possible region for the hidden gas clumpuscules. 



to cover less than ~ 1% of the sky; this means that the maxi- 
mum optical depth should be ~ 0.01 assuming a Galactic halo 
completely made of gaseous clumpuscules. Therefore, we can 
interpret our optical depth limit as the upper limit of the contri- 
bution of turbulent gaseous structures with Rdiff (1. 25pm) < Rd, 
expressed in fraction of the halo (right scale in Fig.[T2lb)). 

Our upper limit does not yet seriously constrain the model 
with clumpuscules, but we can extrapolate these results to de- 
fine a strategy to reach a significant sensitivity. Our monitored 
stellar population is dominated by red giant stars that would give 
a lower scintillation signal than smaller stars from the main se- 
quence. This could only be compensated by achieving an ex- 
cellent photometric precision on the red giants. An easier way 
to improve a hypothetic signal would be to use V passband 
instead of /. In our test, the use of / filter was imposed by 
the choice of the SOFI detector, the only one available with 
a fast readout. A ccording to the ESO exposure time calculator 
dESO-S OFI 2007 ), the same precision we obtained in J for the 
red giants can be reached for AO stars in V with the same ex- 
posure time (10s). Therefore, we can extrapolate that an expo- 
sure of ~ 10 6 star x hour (about 100 times more than our test) 
obtained with the same type of telescope (NTT) using the V 
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passband (around the maximum of stellar emission) should pro- 
vide enough sensitivity to significantly constrain the turbulent 
gas component of the Galactic halo. 

8. Limits on the gas structuration in the nebulae 

For our study of known nebulae, the distance zo to the gas is 
known, but not the star-by-star z\ distances. We first make the 
hypothesis that the stellar population behind the clouds (whose 
light is absorbed and diffracted) is the same as the population 
which is not - or much less - obscured (the so-called control 
population, see Fig. [7]). Toward B68 the red giant stars are distin- 
guished from the main sequence stars with our own / image, also 
taken with the NTT-SOFI detector, through the K s versus (J-K s ) 
diagram. Using Eq. ®, we established the K s to 0^ ax relation, 
providing the maximum angular radius of a main sequence or 
red giant star located beyond 4kpc with K s apparent magni- 
tude (Fig. [TT] right). We conservatively use this angular radius 
in Eq. ([8]) to estimate R™Jff The study of the population behind 
the dust is complicated by the fact that the stars are obscured. 
In the case of B 68 we were able to correct the apparent magni- 
tudes for the absorption in K s band; we used the Ay absorption 
map from ( Alves et al. 2001 ) and the relation A K /A V = 0.089 
( [Ojha et al. 2000[IGlass et al. 19991) to extract A K and deduce the 
dereddened magnitudes K s . Using the relation of Fig. [TT] (right) 
with these corrected magnitudes, we extract 0^ ax for each star 
behind the dust. 

Because there is no absorption map for bcl31 and the 
Circinus nebula, we conservatively used the highest possible 
value of 0f ax ~ 13 x 0(R Q at lOkpc) (corresponding to a B3 
type star at 4.5 kpc). 

Figure [13] shows the cumulative R™jff distributions and the 
upper limits on T2.\6nm(Rd) for the obscured regions of B68 and 
bcl31 (the search regions from Fig. and over the complete 
field of the Circinus nebula, which has indistinct boundaries. 
The best limits are naturally obtained toward B68, as a conse- 
quence of the better knowledge of the stellar sizes. The shape of 
the cumulative R™jff distributions differs from the one obtained 
toward the SMC, because the contribution of (big) red giant stars 
is smaller. 

We can now interpret our limits on Rdiff to constrain the 
global structure of the nebulae, and specifically put limits on the 
existence and the structure of local turbulent dense cores (cells) 
within the nebulae (see lLada et al. 20071 and IRacca 2009b . 

- Probing the global structure? 

If the nebula is a simple object described by a single 
Kolmogorov turbulence law, characterized by an outer scale 
of L z ~ 17000 AU (maximum depth of gas crossed by the 
light in B68), then the minimum expected value of Rdiff is 
deduced from Eq. (Q]), where cr 3n is limited by n max the max- 
imum molecular density. This maximum is estimated to be 
2.61 x I0 5 cm~ 3 for B68 (IHotzel et al.1 2002). Using this con- 
straint and the maximum outer scale value in (Q]) gives for 
B68 



Rdiff (2A6fim) > 
2A6fim 



263 km 



1 /j in 



> 3. x 10 b km. 



11000AU 



10AU 



2.61x10 s cm~ 3 



10 9 < 



,-3 



(13) 



For cbl31, the same calculation based on data from 
dBacmann etal.1 2000) gives R diff (2A6fim) > 4.2 x I0 6 km. 




10 



20 40 60 100 200 400 

R d (km) 




40 60 



R d (km) 



Fig. 13. (a) N*(Rd), the number of directions with no turbulent 
structure of Rdiff(2A6/im) < Rd along the line of sight toward 
the obscured regions of B68 (green), cbl31 (blue) and toward 
the Circinus nebula (red). 

(b) The 95% CL maximum optical depth of structures with 
Rdiff (2A6/im) < R d . 



These large diffusion radii, which are owing to the weakness 
of the power spectrum at small scales, are much too large 
to induce any observable scintillation effect on any type of 
stars. Indeed, because 6 S > 6(R Q at lOkpc), Eq. ^ gives 
R s /R re f > 3100. at A = 2.16yum, and the expected modula- 
tion index is completely negligible (out of scale in Fig. [3]). 
Probing local sub-structures? 

In contrast, local turbulent dense cores with much smaller 
diffusion radii could potentially produce observable scintil- 
lation. Therefore our upper limits on scintillation can be in- 
terpreted as upper limits on the existence of turbulent cells 
with Rdiff(2A6jim) < 350 km within the volume of the neb- 
ula. Using Eq. (Q]), it can be interpreted as limits on the dis- 
tribution of the product 



0~3n 



10 9 



cm 



r][ 



10 AU 



R d iff(2A6fim) 



663 km 



(14) 



along the lines of sight; thanks to the small exponent of L z , 
a rough hypothesis for the core size would allow one to ex- 
tract upper limits on the frequency of cores with a density 
dispersion higher than a given value cr 3 „ . 
No conclusion can be infered for structures with Rdiff > 
350 km, because they cannot produce a detectable scintilla- 
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tion in our sample because of our limited photometric preci- 
sion on the small stars. 



9. Conclusions and perspectives 

The aim of the test was to study the feasibility of a systematic 
search for scintillation. We were lucky enough to find a stochas- 
tic variable light-curve that is compatible with a scintillation ef- 
fect; but considering the low probability of such an event, which 
is related to the low density of the nebulae, a program of syn- 
chroneous multicolor and multi-epoch observations is necessary 
to obtain a convincing signature of the effect. 

From our search for invisible gas toward the SMC, we con- 
clude that an ambitious program using a wide field, fast readout 
camera at the focal plane of a > Am telescope should either dis- 
cover turbulent gas in the halo, or exclude this type of hidden 
baryonic matter. With such a setup, significant results should be 
obtained with an exposure of ~ 10 6 star x hour in V passband. 

The hardware and software techniques required for scintil- 
lation searches are available just now, and a reasonably priced 
dedicated project could be operational within a few years. 
Alternatives under study are the use of the data from the LSST 
project and from the GAIA mission. If a scintillation indication 
is found in the future, one will have to consider a much more am- 
bitious project involving synchronized telescopes, a few thou- 
sand kilometers apart. Such a project would allow one to tem- 
porally and spatially sample an interference pattern, unambigu- 
ously providing the refractive length scale R re f, the speed, and 
the dynamics of the scattering medium. 
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Appendix A: Connection between the diffusion 
radius and the gas structuration 

We consider a gaseous medium characterized by a cell with size 
(L x , L y , L z ). An electromagnetic plane wave propagating along z 
and crossing the medium is distorted. The distorsion is caused 
by the variation of the column density along the propagation of 
the light and can be described by a two-dimensional phase delay 
(p(x,y). We use the phase structure function to characterize the 
variations of the phase as follows: 

D+(x,y) = < (0(x' + x,y>+y) - <P(x' ,y')f > (A.l) 
= 2[< <f> 2 (x',y') > - < 0(Y + x,y' +y)<f>(x',y') >] 
= 2[^(0,0)-^(x,v)], 

where ^(x, y) is the correlation function of the screen phase and 
is related to the phase spectral density S(p(q x ,q y ) (power spec- 
trum per volume unit) through the Fourier transform as 

&(x,y) = J J S^q x ,q y )e 2ni{xq ^ q y ) dq x dq y . (A.2) 

Substituting equation (IA.2I) in equation (IA.1I) we obtain 

Dt(x,y) = 2 J J S^qyXl-e^^^q^qy. (A.3) 

The relation between the two-dimensional spectrum of the 
phase, S^q^qy) and the three-dimensional spectrum of the 
number density fluctuations has been derived for a plasma by 
Lovelace (Lovela ce 19701 lLovelace et al. 19701 ). Here we extend 
the concept to the optical wavelength for a medium of molecular 
gas: 



S^(q x ,q y ) = 2nL z 



A 



(A.4) 



where S^ n (q x , q y , q z ) is the spectrum of the density fluctuations 
of the molecular gas in three dimensions (Sn = n- < n >), 
a is the average polarizability of the molecules and L z is the 
thickness of the medium along the line of sight. 

Assuming the gaseous medium to be isotropically turbulent, 
the 3D spectral density obeys a power law relation within the 
turbulence inertial range: 

Sssniqx, q y , qz) = C 2 n q- p L~ l ut <q< , (A.5) 

where q = ^Jq 2 + q 2 + q 2 , /3 = 11/3 is taken for the 
Kolmogorov turbulence, C 2 is the turbulence strength parameter 
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and L out and L in are the outer and inner scale s of the turbulence 
resp ective ly. By substituting equations (IA.5I) and dA.4b in equa- 
tion (IA.3K we compute the phase structure function in the polar 
coordinate system: 

D^r) = 2Cl(2n)L z (^^-f P° C q*(\ - e 2 ™« cose )d6qdq. 
A Jo Jo 

The integration results in 



D,(r) = 2C 2 n (2nffmJ—\ ^ 



(A.6) 



where 



f(fi) 



-f 

Jo 



s l - p (l-J (s))ds = 



For the Kolmogorov turbulence f(J3 =ll/3)~ 1.118. 

We define the diffusion radius Rmff as the transverse scale 
for which DJ^Rdi/f) = 1 rad. Rdiff is directly obtained from 
equation JA.6t : 



Rdiff = [2C 2 n (2nf +4 f(J])Lz<x 2 *- 2 ] 



1/(2-/3) 



(A.7) 



We will now link the turbulence parameter C 2 to the dispersion 
of the density fluctuation. From Parseval's theorem the disper- 
sion of the density fluctuation is equal to the auto-correlation 
function at origin: 



(A.8) 



On the other hand the auto-correlation is equal to the integration 
over the spectrum in the Fourier space 



m= S 3Sn (q)d 3 q. 



(A.9) 



Using equation (IA.5K the dispersion of the density fluctuations 
is related to C 2 as 



Appendix B: Limitations of the photometric 
precision 

In our photometric optimization, we found that the photomet- 
ric precision on unblended stars was limited by several factors. 
First, the PSF of very bright stars significantly differs from the 
ideal Gaussian and the photometric precision on these stars is se- 
riously affected. Second, the best photometric precision does not 
follow the naive expectation assuming only Poissonian noise. 
For stars fainter than K s = 14.6 (or / = 16.6) we were able to 
reproduce the behavior of the PSF fit x 2 and of the fitted flux 
uncertainties by assuming that the relative uncertainty cry on 
a pixel content results from the combination of the Poissonian 
fluctuations of the number of photoelectrons Nj*, and of a sys- 
tematic uncertainty (i.e. not changing with time, but depending 
on the pixel) due to the flat-fielding procedure: 

l,j L hj J 

where Nj* is the total number of photoelectrons (from the fitted 

star and the - usually dominant - sky background), and is 
the uncertainty of the flat-field coefficient for pixel (i, j). From 
the comparison of two flat-fields taken at different epochs, we es- 
timated that ~ 0.63% (both in Ks and in /). A third source 
of noise in / comes from a residual fringing; we measured that 
fringes contribute to a ~ 40% increase of the Poissonian fluc- 
tuation ^JnJ € j in a domain size comparable to the PSF width 

because of their small spatial scale. 

We conclude that the contribution of the Poissonian noise 
dominates the error when the flux from the star plus the to- 
tal sky background within the PSF fitting domain is small (~ 
2000 ADU = 10600y£). For brighter stars or for stars located in 
a dusty environment (producing a large sky IR background), the 
second - systematic - term dominates. 
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AnC: 



p-3 



mf-'-i?- 3 ). (a. io) 



Assuming L z = L out » L in , for j3 = 1 1/3, we estimate C 2 as 



6nL 



2/3 



(A.ll) 



Finally by substituting equation dA.111) in dA.7lh Rdiff is ob- 
tained in terms of the cloud's parameters and the wavelength 
as 



Rdiff = 231.1 km 



(Xh 2 



1 jjin 



10AU 



0~3n 



10 9 cm~ 3 



(A.12) 



where an 2 = 0.802 x 10 24 cm 3 is the polarisability of molecule 
H 2 . 
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